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Abstract. In many-body and other systems, the physics situation often allows one to 
interpret certain, distinct states by means of a simple picture. In this interpretation, 
the distinct states are not eigenstates of the full Hamiltonian. Hence, there is an 
interaction which makes the distinct states act as doorways into background states 
which are modeled statistically. The crucial quantities are the overlaps between the 
eigenstates of the full Hamiltonian and the doorway states, that is, the coupling 
coefficients occuring in the expansion of true eigenstates in the simple model basis. 
Recently, the distribution of the maximum coupling coefficients was introduced as 
a new, highly sensitive statistical observable. In the particularly important regime 
of weak interactions, this distribution is very well approximated by the fidelity 
distribution, defined as the distribution of the overlap between the doorway states 
with interaction and without interaction. Using a random matrix model, we calculate 
the latter distribution exactly for regular and chaotic background states in the cases 
of preserved and fully broken time-reversal invariance. We also perform numerical 
simulations and find excellent agreement with our analytical results. 
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1. Introduction 

In open quantum systems, strength function phenomena pQ give structural information 
about the system itself and about the excitation mechanism. Here, we address statistical 
features of the doorway mechanism which can be defined as follows: there are one or 
several somehow "distinct" and "simple" excitations whose amplitudes are spread over 
many "complicated" states. In a many-body system, collective excitations are often 
distinct, because all or large groups of particles move in a coherent fashion. As compared 
to the complexity of the other, non-collective excitations, these states can be interpreted 
in the framework of a ted in the framework of a simple, typically semiclassical, picture. 
The distinct states act as "doorways" to the background of the complicated states [JJ[2]. 
Mostly, the statistical features of the latter are chaotic. The strength function has Breit- 
Wigner shape, largely independent of the statistics of the background states. The width 
characterizing the Breit-Wigner strength function is referred to as spreading width [JJ. 

The doorway mechanism is found in a rich variety of systems, comprising atoms and 
molecules [3], as well as atomic clusters, quantum dots and, more generally, mesoscopic 
systems [H [5j [6]. Nuclear physics provides particularly beautiful and well-studied 
examples, such as Isobaric Analog States and multipole Giant Resonances [Tl 171 151 l9l [TO] . 

What is a suitable theoretical interpretation of the Breit-Wigner shape? 
Although the simple picture for the distinct excitations captures the main physics, 
it is important to realize that these states are not eigenstates of the real quantum 
Hamiltonian. Similarly, the statistical models for the background states do not describe 
eigenstates either. Thus, if we use the simple picture for the distinct states and the 
statistical model for the background states as a basis of the Hilbert space, there must 
be a non-vanishing interacting between these two classes of states. Rediagonaliziation 
then yields proper eigenstates of the model Hamiltonian. Averaging over the background 
states, one obtains the local density of states around the energy of the distinct state 
or states in the simple picture. The local density of states is once more of Lorentzian 
or Breit-Wigner shape with a spreading width that is — depending on the particular 
situation — closely related to or identical with the above mentioned spreading width 
in the strength function. It can be viewed as a measure for the quality of the simple 
picture describing the distinct states: the smaller the spreading width, the closer is this 
picture to the physics reality. 

The strength of the interacting between the two classes of states uniquely determines 
the spreading width and, equivalently, the size of the overlap between the distinct state 
in the simple picture and the true eigenstates of the model Hamiltonian. These overlaps 
are of course the coupling coefficients when expanding the true eigenstates in the above 
mentioned basis. Recently, a new statistical observable was introduced: the distribution 
of the maximum coupling coefficients [UJ. The first two moments of this distribution 
were already studied in Ref. [12], but with assumptions not valid in our context. In 
Ref [IT] , however, the full distribution is addressed. Importantly, its shape sensitively 
depends on the interaction strength. Moreover, it is an especially well-tailored measure 
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to investigate weak interactions. 

Here, we present exact results for the distribution of the coupling coefficients to a 
distinct state in the framework of a random matrix model. In the particularly interesting 
regime of weak interactions, this distribution coincides with the distribution of the 
maximum coupling coefficients. 

The article is organized as follows. After properly posing the problem in Sec. [21 
we calculate the distribution exactly for regular and chaotic background, respectively 
in Sees. [3] and HI We discuss our results in Sec. El 

2. Posing the Problem 

In Sec. 12.11 we present the random matrix model for the doorway mechanism. We 
introduce and define the distribution of the maximum coupling coefficient in Sec. 12.21 

2.1. Doorway Mechanism in a Random Matrix Model 

The model to be discussed here stems from nuclear physics |T] and is also often used 
in other fields pL3j . For the convenience of the reader and to define our notation, we 
compile its salient features. As we are aiming at a random matrix model, it is convenient 
to choose from the beginning a proper basis of the full Hilbert space such that we can 
represent the Hilbert space operators by matrices. Introducing a cutoff, their dimension 
is finite. Eventually this cutoff effect is removed by taking the matrix dimension to 
infinity. We nevertheless use the Dirac notation for the wave functions, even though 
they are finite-dimensional vectors. 

The total Hamiltonian H consists of three parts, the Hamiltonian H s for the K 
distinct states which become the doorway states, the Hamiltonian H b describing the N 
background states, where N will eventually be taken to infinity, and the interaction V 
coupling the two classes of states. Hence, we have 

H = H s + H b + V 

K N K N 

= j2 E sj\s j )(s j \+J2 E ^^ ( y ^i s i><M + h - c ) • w 

j=l u=l j=l v=\ 

For the matrix elements of the interaction, we make the assumptions (sj\V\sk) = 
(K\V\bfi) = and (b u \V\sj) = V V j for any j, k, /x, v. Often, there is only one relevant 
doorway state or the spacing between the doorway states is much larger than their 
spreading widths. We focus on these cases and consider only one doorway state by 
setting K — 1, \si) — \s) and Vj U = V u . 

The eigenequations for the uncoupled Hamiltonians are 

H s \s} = E s \s) and H b \b v ) = E u \b u ) (2) 

Due to the interaction V the doorway state is not an eigenstate of the Hamiltonian H. 
We denote the eigenstates of the full Hamiltonian H by \n). The eigenequation to be 
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solved is 

H\n) = E n \n) (3) 

Resembling the situation in most systems, we put the doorway state \s) in the center of 
the background spectrum. It interacts with the surrounding N states. Without loss of 
generality, we may set E s = 0. 

The exact eigenstate of H which evolves from the doorway state in the presence of 
the interaction is referred to as |0). We expand the n— th eigenstate of H in the basis 
spanned by \b v ) and \s) as 

N 

\n) = c ns \s) + y^ y c nv \b v ) (4) 

where the coupling coefficient c ns is the overlap between the doorway state \s) in the 
non-exact picture for this distinct state and the n-th exact eigenstate |n) of the full 
Hamiltonian. We are interested in the statistical features of these coupling coefficients. 

We have to solve our model for c ns . The action of the full Hamiltonian H on the 
eigenstate |n) yields on the one hand 

tJV \ N 

E s C ns + KCnv \\s) + J2 {*» V » + E » C ™) l 6 "> ' ( 5 ) 
u=l J v=\ 

On the other hand we have 

TV 

H\n) = E n c ns | s) + E n c nv \ b v ) . (6) 
Equating these two expressions, we find 

c nu — j-, j-, c ns j {' ) 

hi n — hi v 

such that 

|n) = c ns (\s) + ^^^Jb„)j . (8) 
Using the normalization of \n), we eventually arrive at 

N n/i 2 V 1 

which is the desired expression for c ns in terms of the matrix elements of H. 

Formula ()9]) holds for < n < N . This expression is still exact. However, since 
the eigenvalues of the full Hamiltonian E n depend on the coupling coefficents V u , it is a 
complicated implicit expression. We proceed further by expanding the exact eigenvalues 
perturbatively in V, 



N 



E n = E u{n) + ^ + J2 lVvi ts +••• (l<n<iV) 

N IT/12 N N |T/|2|T/|2 

^^-Ef-EE^f. do) 



|2 
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where the eigenstate |n) of the full Hamiltonian to eigenvalue E n has evolved from the 
eigenstate \b u (n)) of the unperturbed Hamiltonian by adiabatically switching on the 
perturbation. 

We now obtain a crucial simplification by setting E n = E v ( n ), i. e. we keep only the 
leading order term in the perturbative expansion of Eq. ( fTOl) . To this approximation 
of |c ns | 2 for ti / the sum in Eq. (Q is completely dominated by the term v = u(n), 
which actually diverges such that |c ns | 2 ~ to first order for all The only overlap 

integral which remains finite is the overlap of the doorway state with itself |co s | 2 . Here 
we set E E s = and no divergence occurs. Using the approximation E n = E v ( n ) 
we have therefore essentially singled out |co s | 2 as the only non- vanishing - and thus 
inevitably maximum - overlap integral of the perturbed eigenstates with the doorway 
state. 



2.2. Distribution of the Maximum Coupling Coefficient 

The new statistical observable introduced in Ref. [11] is the distribution of the maximum 

Cmax = max(|c ns |,0 < n < N) (11) 

of the overlaps between the eigenstates of the full Hamiltonian and the distinct state, 
that is the doorway state \s). In order to obtain it, we have to average in a suitable way 
over the interaction matrix elements and over the Hamiltonian modeling the background 
states. For the time being it suffices to denote this average by square brackets. Later 
on we give a precise definition. Hence, the distribution in question is given by 

Pmax(c) = (S(C- C max )) . (12) 

On the other hand, the distribution of overlap between the evolved doorway state and 
the unperturbed doorway state reads 

po(c) = (5(c-\c 0s \)) . (13) 

Setting E n = E v ( n ) amounts essentially to the approximation 

Pmax(c) « p (c) . (14) 

This approximation is certainly good for small interactions or, more precisely, as long 
as the mean coupling strength is an order of magnitude smaller than the mean level 
density of the background states. Our numerical simulations will strongly corroborate 
this statement. Hence, we focus on Po(c), which can be treated analytically. 

The statistics of the interaction matrix elements can only have minor impact on 
the resulting distribution. Hence, if not stated otherwise assume that the interaction 
matrix elements are Gaussian distributed random variables. We have to distinguish two 
cases. The total Hamiltonian H can be time-reversal non-invariant or time-reversal 
invariant, where we disregard spin degrees of freedom. In the first case, labeled by the 
Dyson index (3 = 2, the interaction matrix elements V u are complex variables, in the 
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second case, labeled (3 = 1, they are real. Introducing the ^-component vector V, the 
corresponding distribution is 

In addition to the behavior under time-reversal invariance, the statistical properties of 
the Hamiltonian Hj,, however, must strongly affect the distribution po{c). Hence, we 
do not specify it yet. As is well known from Random Matrix Theory, the parameter 
governing the physics is 

a = jM « (16) 

where D is the mean level spacing of the background states in the center of the 
band [HI3]. The distribution Pi(V) is chosen such that A is independent of 0. 

Technically, it is more convenient to work out the probability density Q(u) of the 
random variable 

1 - II/ I 2 

The relation between the two distributions reads 

du 



Po(c) = Q(u) 



dc 



. 3 

u=l/c' 



(18) 

u=l/c 2 



Thus, once Q(u) is known, po{c) follows immediately. 

We now use our statistical assumption that the interaction matrix elements V u are 
Gaussian distributed. We write the distribution Q(u) in the form 



Q( U ) = divmv) u u-i-r^ , (19) 




where d[V] is the product of the differentials of all independent variables in V. The 
square brackets with index N denote an average over the iV background states, that is, 
over the Hamiltonian Hf,. For the calculation of the averages, it is helpful to write the 
distribution Q(u) as the Fourier transform 

+oo 

Q(u) = — dkexp{ik(u-l))R{k) (20) 



2tt „ 

— oo 

with the characteristic function 

R(k) = [ d[V]P t {V) ( exp ( -ik > ] ) ) . (21) 

N 

After rescaling V u = y v /\E v \ we obtain the alternative expression 

-0N 




R(k) = y J d[y] exp {-iky ] y) 

detH b fe^(-^H^yjJ , (22) 
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where the vector y has real entries for (3 = 1 and complex ones for (3 = 2, respectively. 
The infinitesimal volume element d[y] is a product of the differentials of all independent 
entries of the vector y. As a generating function R is normalized to R{0) = 1. 

In the sequel, we calculate the expressions (|22|) for generic choices of the 
Hamiltonian H b governing the dynamics of the background states. 



3. Regular Background 

The doorway state is embedded into a regular background, if the eigenvalues E v of H b 
do not repel each other. The distribution of the background Hamiltonian then factorizes 
according to 

N 

P b (H b ) = Y[p b {E v ) . (23) 

In order to keep the discussion most general we use for the interaction matrix elements 
a general factorizing distribution 

TV 

Pi{V) = l[pi(V v ) , (24) 

instead of the Gaussian distribution introduced before Eq. ffT5]) . We keep the reasonable, 
physically motivated assumption of statistical independence of the interaction matrix 
elements but relax the global orthogonal ((3 = 1) or unitary ((3 = 2) invariance 
of the interaction matrix elements, implicit in the measure Eq. (Tl5|) . For complex 
coupling matrix elements we assume in addition that the distribution p^ is U(l) invariant 
Pi(V) = pi(\V\). We asign to complex coupling matrix elements with this invariance the 
Dyson index (3 = 2 and to real coupling matrix elements the Dyson index (3 = 1. 

A straightforward calculation reveals that the characteristic function (I2T]) factorizes 
as well and becomes an iV-th power of a single integral, 

R(k) = e Nlnr{k) (25) 

+oo 

r(k) = J d%\ Vl {z) J dE Pb (E)exp (--^\ • (26) 

— oo 

As we are interested in the local scale set by the mean level spacing D of the background 
states, the distribution po(c) should not be sensitive to the particular choice of the 
distribution p b , as long as it does not contain scales competing with the mean level 
spacing D. The simplest choice is 

. , , \E\<VN/2 

*W = < n \E\>VN/2, (27) 




where D = 1/ y N and y N = ND is the length of the background spectrum. The 
following calculation is similar to the one described in the Appendix B of Ref . [H] . We 
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perform the integral over the background distribution in Eq. 

-L-^O OO 

f-ik\z\ 2 \ 2 f exp (-ik\z\ 2 u 2 ) 

dEp b (E) exp — — — = -= / du- 



E 2 J v/jV J u 2 

2/VN 



n^oo 2 f du ex V (-ik\z\ 2 u 2 ) 



N J u 





2 



= 1 + 2|^. (28) 
In the last equation we used an integral identity of the Fresnel type 

oo 

id f ^exp(- W) = to 

- ■ £ \l A r\i 



da J u 

o 

We find for the characteristic function 

R(k) = exp (^VinkNm^j (30) 

m x = [ d p [z] Pi (z)\z\ . (31) 



We observe that the distribution of the interaction matrix elements Pi enters only via the 
expectation value mi as defined in Eq. (130]) and not via the second moment m 2 = v 2 . As 
pointed out after Eq. (JT6l) pi is chosen such that the mean coupling strength, as defined 
through v, is independent of j3. This means that mi in general is different for real and 
for complex coupling. We write mi = a@v, where now ap depends on Dyson's index and 
on the distribution p^. For instance for the Gaussian distribution we find 



(G, = ^ « 0.80, = 1 
1 % w 0.89 , (5 = 2 . 



Using the definition of A in Eq. f|T6|) we finally find 



R{k) = expl^afiXV^Kk) . (33) 



The reader can easily convince herself /himself that other reasonable choices for Pb{E) 
such as a Gaussian distribution yield the same functional form as in Eq. fl33|) . The 
Fourier transform fl20l) results in 



Q(u) = (^^^(-(a.A) 2 -^) . (34) 
Using the relation f|T8|) . we eventually arrive at 

Po( c ) = 77— ex P (~ (M) 2 T— T ) ' ( 35 ) 



As anticipated, the interaction strenght enters the distribution only via the dimensionless 
ratio A = v/D. The distribution pi enters via the the factor ap defined in Eq. f l30|) . 
It is interesting to see that ap does not only depend on the distribution but also on 



Exact Coupling Coefficient Distribution in the Doorway Mechanism 



9 



Po (c) 




0.0 0.2 0.4 0.6 0.8 1.0 

Figure 1. Plot of the the distribution function po(c) for a regular background for 
real interaction matrix elements (dotted) and for complex matrix element (full line) 
for three different values of the mean coupling strength A = 0.1 (blue), 0.5 (green), 2 
(red) . The difference between real and complex coupling increases for strong coupling 
A. 

the symmetry factor 0. This means that for a regular background and for constant 
interaction strength A the distribution p distinguishes between real interaction and 
complex interaction. As we will see in the following this does not happen for a chaotic 
background. Therefore here opens - at least theoretically - the possiblity to distinguish 
between a regular and a chaotic background dynamics through the doorway state. 

Assume we can experimentally manipulate the interaction between the Doorway 
state and the background such the interaction matrix elements change from real to 
complex, for instance by switching on a magnetic field. For fixed mean interaction 
strength and for a chaotic background dynamics the distribution p will be invariant, 
whereas for a regular background it will change. 

In Fig. [Upo is plotted for three different values of the mean coupling strength. We 
see that due to the numerical similarity of a\ and a 2 the plots for real and for complex 
couplings are almost the same. The fact that for a Gaussian distribution a\ pa a 2 seems 
to be rather accidental. However for p^ being a semicircle (SC) distribution we find 
a^ SC ^ pa 0.85 and c4 SC ' ) ~ 0.93 and for a uniform (U) distribution we find = a/3/2 
pa 0.87 and a 2 — a/8/3 ~ 0.94. In other words, for physically reasonable distribution 
functions pi we have a\ pa a 2 , which encumbers the experimental possibility sketched 
above. 

4. Chaotic Background 

The dynamics of the background states is usually chaotic. The N x N Hamilton 
matrix modeling the background states has then to be chosen from a Gaussian 
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random matrix ensemble. Again, we have to distinguish time-reversal invariant and 
non-invariant systems, that is, the cases labeled by the Dyson parameters (3 = 1 and 
(3 = 2, respectively. The Hamiltonian H b is from the Gaussian orthogonal ensemble 
(GOE) for (3 = 1 and from the Gaussian unitary ensemble (GUE) for (3 = 2 with 
variance w 2 of the diagonal elements 

P b {H b ) ~ expf-^tr# 6 2 J (36) 

As already said, y is an iV-component vector which has real or complex entries for (3 = 1 
and (3 = 2, respectively. 

In Sec. 14.11 we reformulate the problem in terms of matrix invariants. This enables 
us to introduce a handy supermatrix model in Sec. 14.21 which we then solve in Sec. 14.31 



4-1- Reformulation in a Rotation-Invariant Form 

At first sight, the best way of tackling the problem seems to start from Eq. (jT9l) and to 
introduce V u = y v /\E v \. The y integration is then simply an integration over a JV — 1 
sphere in a real or complex space. This leads to 

QM = %feil^/ |det ^ e -»-A.., (37) 



2 {2vH/(3)-P N / 2 V 01 /n 
where the N x N matrix 

G = diag((u- l)/3/2t» 2 , 0, . . . , 0) (38) 

has rank one. Furthermore, 

a « - wm (39) 

is the volume of the above mentioned sphere. Unfortunately, the remaining ensemble 
average is only feasible for the GUE, we carry it out in Appendix A For the GOE, the 
calculation is hampered by the modulus of the determinant. 

To adress the GOE case and to have a method that is capable of handling both cases, 
GOE and GUE, in a unifying way, it turns out necessary to cast the ensemble average 
into an invariant form. To this end, we start from Eq. (|22|) . perform the integration 
over the vector y and obtain 

R(k) = (det(H 2 + 2tv 2 k/(3)-^ 2 \ det H b \ p ) N . (40) 

Now we make the observation, that the N x iV-matrix average can be expressed as an 
(N + 1) x (JV + l)-matrix average as follows 



(det(ff 2 + 2iv 2 k/(3)-P /2 \ detH b f) N = lim L Nf} w^ N V2nw 2 

e— >0 

(e + 2iv 2 k/(3f 12 Ui5(H)det(H 2 + 2iv 2 k/(3Y^ 2 ) , (41) 

\ / N+l 

where H is a (N + 1) x (JV + 1) matrix and the ensemble average is over an 
(JV + 1) x (JV + l)-matrix GOE (GUE) ensemble. The derivation of Eq. flUD, which is 
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crucial for the calculation, is sketched in App. Appendix B On the right hand side, the 
inconvenient modulus of the determinant has disappeared. The combinatorial factor 
L NI3 is given by 

2r(l + (iV + l)/2) 

L Nfi = { vMiV + 1) ' P . (42) 
M , [3 = 2 



We express the determinant on the right hand side of Eq. (1411) as a Gaussian integral 
over a real ((3 = 1) or complex (/3 = 2) N + 1 vector y. 

det(H 2 + 2iv 2 k/(3 + e)- pl2 = J ^ exp {-y\H 2 + 2iv 2 k/(3 + e)y) (43) 

and plug subsequently Eqs. (HOI) . ( ]4~Ti) and ( )43l) into Eq. ( 120]) and write the integral over 
the iV + 1 vector y in radial coordinates. We find 

L Np n N+l ( f3w 2 x 

Q(w) - 



27^/ 2 V 2V 2 71 



oo 

lim J dxx^ N+1 ^ 2 - l F N+l (l + ^) 





oo 



dk(e + tk)^ clk(u _ x _ x) ^ 
2n 

where we introduced the function 

F N+1 (g) = V2^(ti5(H)e-^ tl5n2 ) (45) 

\ / N+l 

with the {N + 1) x (N + 1) matrix G = diag (g - 1, 0, . . . , 0). 
After combining the various constants we arrive at 

'(3w 2 V N/2 d?' 2 



Q(u) 



r(/?/2) \ 2v 2 J diu-iy/ 2 

■2 

2~ 



dxx ^)/ 2 ^F N+l ( 1 + ^f] 5(u-l- x 



V 

where we introduced formally the fractional derivative 



oo 



dx?l 2 v ' e->o2vr 



We can evaluate the fractional distribution 



(46) 



' r 5(x) = lim -!- / rfA;(e + ikf /2 e ikx . (47) 



2 d r r e+^r(f + l) 
— — Im hm v % /9 ' 48 
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for arbitrary /3. For /3 — 1, 2 we obtain 

d 2r(3/2) 



d^ 2 
dx?l- 



;S(x) 



dx 7T^/x 

l 5 ^ 



e(x) (3 = i 

= 2. 



(49) 



Here we see that the GOE case is more complicated than the GUE case. 

For the GUE the fractional derivative disappears and we find without further 
problems 



d ( f3(u-l)w 2 \ N ~ ( P(u-l)w 



(50) 



Q{U) du\ 2v* 
For the GOE we obtain an integral expression 

o 

The remaining task is in both cases (GOE and GUE) the calculation of F^ + i(g). 

4-2. Mapping onto a Supermatrix Model 
Using 

1 1 

tr 6(H) = -Im tr 

7T H — it 



1 _ d det(H + j) 
— 1m ■ 



2tt dj det(H - j - ie) 



(52) 



3=0 



the ensemble average F^+iig) defined in Eq. (145]) can be expressed via standard 
techniques supersymmetric matrix integral 

1 d i / o \W 2 



7V+1 



is) 




Im 



- 



2vr dj y/g \ 1 + g 
d[rj exp 



i>gw z 



St, ( r + i±£j) ] Sdet~ 1 (r- + J) 



1-5 



exp 



1 



2w 2 



Stra 2 Sdet-^/ 2 (a- + r^ 



(53) 



3=0 



Here a and r are 2x2 (GUE) respectively 4x4 (GOE) supermatrices of the form 

, GUE 



CLi A* 

Ai ia 2 

( ai a 2 A* — Ai \ 

a 2 ^3 Ag — A2 

Ai A2 ici4 

V A* A2 ia 4 J 



GOE . 



(54) 
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The matrix entries in latin letters denote real commuting integration variables. The 
matrix entries in greek letters denote complex anticommuting integration variables. 
The infinitesimal volume elements d[r\ and d[o~] are products of the differentials of 
all independent integration variables. The integration domain of the real commuting 
variables is the real axis. The matrix J is a 2 x 2 (GUE) or a 4 x 4 (GOE) diagonal 
supermatrix with entries J = diag (j, — j) (GUE) and J = diag (j, j, — j, — j) (GOE). 
Due to the broken rotation invariance of the original matrix model (JT|) the resulting 
supersymmetric representation ( 1531) is a two-matrix model. 

We wish to evaluate the a-integral by a saddle-point approximation and to 
calculate the r integral exactly afterwards. It is well known [15] that the a integral 
Kn{t) yields for large N in the saddle-point approximation 



K n (t) = J d[o-]exp (-^ Strfj2 ) Sdet-^V" 



2 iflnStrr \ ^ ( 1 



exp |__ St ^ + ^_ +0 -| , (55) 



where D = \J (3it 2 w 2 / (2N) is the mean level spacing in the center of the band. However 
this approximation is only valid, if Str r itself is of order of the mean level spacing. Since 
the integration domain of r is the whole real axis, this is not automatically guaranteed. 
A necessary condition is that the variance of the Gaussian in the second line of Eq. (1531) 
is itself of order of the mean level spacing, i. e. the r integral in Eq. ( )53l) is essentially 
localised to a small window of width D around zero. Consequently we must require 

~ N + 0(1) . (56) 

And therefore g should scale as iV for large N. The dimensionless coupling strength is 
given by A = v/D. Since u it is of order one, we obtain for g in the GUE case 

2(u-l) 

Fortunately this is exactly the scaling behaviour we need to apply the saddle-point 
approximation. In the GOE case Eq. ( IBTl) holds in any Intervall u c < x < (u — 1), 
where u c is an infrared cutoff in the integral, which is small compared to one but large 
compared to the mean level spacing, i.e. in the limit N — > oo in the whole integration 
domain of the x-integral in Eq. flBTl) . In conclusion we can apply the approximation 
(|55|) both in the GUE as well as in the GOE case. 

4-3. Remaining Matrix Integration and Final Result 



Plugging Eq. (155j) into Eq. (1531) we obtain after a simple shift 

~ . . ft? d 1 / 2 Y 1 "' 2 (2m 
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, (58) 



j=0 



where we also employed that g ~ N 1. Now the derivative with respect to the source 
term can be performed 

~ , , 1 / o N W 



1 + Ry^l m ( tr 



where we used the identity 

1 

tr 



H + z 



N 



l_d 

2dj 



c/[r ] e -^St r ( r+ ,-J) 2 sdet - /3 7V/2 r - 



(59) 



(60) 



j=0 



for a complex number z with negative imaginary part. The remaining average can be 
calculated by employing techniques of standard analysis. We use 



Im — / dx 

2tt 



exp (—x 2 /2) 



x — tr 



sgn(r)W — exp I — erfc \ — 




which holds for r£l. We obtain 

1 

Im ( tr 



TJ i BlTTW 2 



7T 

2«7 



exp 



(3N 



erfc 



' f3N 



(61) 



(62) 



Finally we obtain 

F N +i(g) 



'27TW 2 



gD 2 \l + g 



0N/2 



1 + — \l "7TTT exp 



/3iV 



erfc 



'/3iV 



2\j f3N 

This result simplifies considerable when we take into account the scaling behavior (J5 
of g. In the large N limit we can write 



(63) 



, (5xw' 
hm — 

7V->oo V 2v 2 



2\ 0^/2 



N+l 



1 + 



(3xw* 



'2ttA 2 



exp 



2x 



+ erfc 



2x 



(64) 



This can now be plugged into Eq. fj46|) to obtain expressions for Q(u) on the scale of 
the mean level spacing. For the GUE we find straightforwardly 
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Figure 2. Figure showing the plots of the analytical results for the GOE (blue curves) 
and for the GUE (red curves) for three different coupling strength (A = v/D) A = 0.1 
(thick curves), A = 0.5 (dashed curves) and A = 2 (dotted curves) 



For the GOE we are left with an integral expression for Q(u) 



Q{u) 



7T 3 A 6 



dx 



2{u-lfJ yi^r 2 
o 



cxp 



7T 2 A 2 



2{u - l)x 



(66) 



The integral can be evaluated further and be expressed in terms of standard special 
functions. Finally we arrive at 



Q(u) 



Po{c) 




7T 2 X 2 C 2 



4(1 -c 2 ) J ' " x \4{l-c 2 ) )\ ' ^ 
where K n is the modified Bessel function of the second kind of order n. In the Fig. (T2J) 
the distributions po(c) of Eq. (IBT|) for the GOE (blue curves) and of Eq. ([55]) for the GUE 
(red curves) are plotted for the values A = 0.1, 0.5, 2 of the mean coupling strength 
A. We see that for small A there is only a minor difference between GUE and GOE 
background. 



4-4- Comparison 

In Fig. [3] the distribution function of the overlap integral |(0|s)| of the evolved doorway 
state with the unperturbed doorway state is plotted for four different coupling strengths 
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Figure 3. Four figures showing the analytical results 1) for a complex coupling of 
the doorway state to a GUE background (full red line) 2) for a real coupling to a 
GOE background (full green line), 3) for a complex coupling to a regular background 
(dashed blue line) and 4) for a real coupling to a regular background (full blue line) for 
four different coupling strength A = v/D: A = 0.05 (upper left panel), A = 0.1 (upper 
right panel), A = 0.5 (lower left panel) and A = 2 (lower right panel). 

A = 0.05, 0.1, 0.5 and 2 and for all types of couplings and background complexities 
considered. These are: 1) complex coupling to a regular background (dashed blue 
line), 2) real coupling to a regular background (full blue line), 3) complex coupling 
to a GUE background (full red line) and 4) real coupling to a GOE background (full 
green line). As a general trend mixing with the background is strongest for a complex 
coupling to a GUE background and weakest for real coupling to Poissonian background. 
However the difference of the distributions for different background complexities is rather 
small. This suggests a certain degree of universality of the curves. One might choose 
other ensembles for the background Hamiltonian, as for instance semi-Poisson|16] or 
transition ensembles. However we expect that for these ensembles which lie between the 
two extreme cases, GUE and Poissonian, their corresponding distributions will also lie 
in the channel between the full red line (GUE) and the full blue line (Poissonian with 
real coupling). For the most interesting case of small A this channel is small. 

On the other hand the distributions are highly sensitive with respect to a change 
in the coupling strength A. 

In Fig. H] we compare the curves for po(c = |(0|s)|) obtained from the analytical 
results (in this case from Eq. ( 167|) ) with Monte Carlo simulations in the case of a real 
coupling to a GOE background. The figure shows p (c = |(0|s)|) for three values of the 
coupling strength A = 0.1, 0.5 and 2. We see fairly good agreement for all three values, 
even for the strong coupling value A = 2. This shows that the approximation implied 
in Eq. ( flOl) and thereafter, is justified far beyond the perturbative regime. 
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Figure 4. Comparison of the analytical curves (dashed lines) obtained from Eq. (p7|) 
with Monte Carlo simulations (full lines) for A = 0.1 (red), A = 0.5 (green) and A = 2 
(blue) for a GOE background and real coupling coefficients. 

5. Discussion 

The distribution of the maximum coupling coefficients in the doorway mechanism has 
been introduced as a new statistical observable. These coupling coefficients, that is, 
the overlaps between the eigenstates of the full Hamiltonian and the doorway state 
are not always at the disposal. However, in situations where they are accessible, 
this distribution provides a highly sensitive measure for the interaction strength. Of 
particular interest is the regime of weak interactions. In this regime, the distribution of 
the maximum coupling coefficients is very well approximated by the distribution of the 
overlap between the evolved doorway state and the unperturbed doorway state. While 
calculating the former seems unfeasible at present, we calculated the latter exactly for 
regular and chaotic background states in the cases of preserved and fully broken time- 
reversal invariance. We performed our calculations in the framework of Random Matrix 
Theory which is well-known to provide reliable models for regular and chaotic systems. 
We also carried out numerical simulations which fully confirm our analytical results. 

Our exact calculations are of general interest for matrix models. We managed 
to reformulate a problem with breaking of rotation invariance in the space of N x N 
random matrices in terms of a rotation invariant problem involving (N + 1) x (N + 1) 
random matrices. This made it possible to map the matrix model in ordinary space 
onto a matrix model in superspace, which we solved by a saddle point approximation 
in the limit of infinite level number. Remarkably, the supermatrix model is in the class 
of two-matrix models which show up in a large variety of situations. 
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Appendix A. GUE Background with Finite Level Number 

We define the average 

(...) 9 n = A N J d[H](...)e'^ tTG ' H2 , (A.l) 

where the diagonal NxN matrix G is defined as G = diag (g', 1, . . . , 1). Her g' is related 
to the parameter of the main text as g' — 1 + 2(u — l)w 2 /v 2 . The integration is over the 
set of all Hermitean NxN matrices, i. e. over the GUE ensemble. The normalisation 
An is chosen such that (1)ijv = 1- The task is to calculate the average 

F N (g') = (detH 2 ) g , N . (A.2) 

A Laplace expansion of the determinant yields 

N 

F N {g')= (- 1 ) sgnM+ssn(a, ' ) (II // -H // -'H")^' ( A - 3 ) 

u,u'£Sjv 71=1 

where Sn is the permutation group. Obviously only terms uj — u/ contribute. 

N 

Fn(9') = £<ni ff Mn)lVr, (A.4) 

weSjv n=i 

It is useful to expand the remaining sum in cycles involving the index 1. For indices 
k n > 1 and k n ^ k m we define the cycles Cn^ as 

C^\ki, . . . k n ) = \Hi kl \ 2 \H klk2 \ 2 . . . \H kn _ lkn \ 2 \H kn i\ 2 

C« = \H n \ 2 (A.5) 

Since the indices 1, ki, . . . k n do not appear in the remainder of the product, we can 
integrate over the remainder separately. This yields 

N-l 

F n(9')= E E (Ci^---A)WWi(l) (A.6) 

The average over the cycles is simple as well. Only terms involving the index 1 yield a 
factor different from w 2 . We ontain 



Exact Coupling Coefficient Distribution in the Doorway Mechanism 19 



N+l 



2w 4 



2 \ N+1 w 2( - n+1 ^ 

(C«(* 1 ,...* B )>= (—J — n>l (A.7) 



The averages are independent of the indices k n . The sum over the indices yields the 
combinatorial factor (N — l)!/(iV — 1 — n)\. Altogether we obtain 

2 



2(JV-1K( — ] F*_ a (l) + 



2 Vy^ +1 ) (iV-1)! 

Evaluating this equation for g' = 1 allows us to replace the sum. After some further 
simple manipulations we finally obtain 



7^1 E^T^SfW.d) . (A.8) 



which is almost our final result. The remaining task is to evaluate the g' independent 
constant Fn(1). This is facilitated by the observation that 

F N (1) = ^(/3w 2 ) N MK N+1 (0,0) (A.10) 

where K N+1 (x,y) = J^n^^ni^^iv) is the standard GUE kernel as defined in 
Eq. (6.2.10) of Metha's book (third edition). The 

<f> n [x) = (2 n n\^r 1/2 ex P (-x 2 /2)H n (x) (All) 

are oscillator wave functions and H n is the n-th Hermite polynomial. The constant 
Kn(0, 0) can also be evaluated 

k ro m 1 m / K iV - 1 )/ 2 )r 2 at odd 

^ ( °' 0) = ^2^\ [(iV-2)/2)!]- 2 iVeven. (A ' 12) 



Of course Kn(0,0) is but the inverse level spacing at the center of the semicircle. 
Therefore limAr^oo Kn(0, 0) / \/2N — 1/tt. We obtain our final result 

N + 1 i ^2 



where 



F N (g>) = (2wYN\K N+1 (0,0)^(^) (^^ + C »)^ 

f 1 N ° dd ( K 1 A\ 

Cn = { 1 + 1/N N even (A ' M) 

The even-odd difference disappears in the large N limit. In the following we set = 1. 
Using Eq. (157)) we find for Q(w) 

2\ at 

w A - AT-l 



Q(u) = ( — ) IC^+iCO, 0)(w— 1 



F ^ 2 ^ + Vte'-^ +JV] , (A , 5) 
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This exact result can be compared with Eq. (150 j) in order to find a differential equation 
for F N+1 (g') 

(^ + (^-i)| ? )f^ = ^(o > o)^(^)" +1 

This differential equation can easily be solved. However the solution 

/ o n N 

71 I I 



N) . (A.16) 



" '.t^tJ y ^(x+i)^y (A - ir) 

is highly complicated. It discourages any attempt to calculate the matrix integral 
Eq. ( 15 3 p for finite N in the GOE case. To make contact with the results obtained 
in the main text, we introduce the function 



Pis^ = \l-n—^2- (A.18) 



2Ng > 

With a change of variables y = p(x) in the integral we can write 



N 

'■V-H = J J, i^rj) + i .(0.0) 



which coincides with Eq. (163 p in the large N limit and for /3 = 2. 



(A.19) 



Appendix B. Derivation of Eq. ( 14 lj) 

We define the function 

G N (z) = (det(iJ 2 + ^)~ /3 / 2 | det l^)^ (B.l) 

where H is a A x A GOE or GUE random matrix and the brackets denote the 
corresponding GOE or GUE average. Gn (z) is an analytic function in the cut complex 
plane C \ M_ . In this Appendix we prove the identity 

G N (z) = L N ^' 2 (ti 8(H) det(H 2 + z)'^ 2 ) , (B.2) 

\ / N+l 

where if is a (N+l) x (N + l) GOE or GUE random matrix. The constant L^p is given 
in Eq. (|42|) . We write the rhs of Eq. (1B.2[) in angle eigenvalue coordinates H — > U~ 1 EU ', 



where £isaiV + lxiV + l diagonal matrix of the eigenvalues E!j of if. Since the 
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average is over an invariant function the integral over the diagonalizing group is trivial. 
The average on the rhs. can now be written as 

N+l N+l 



z] -fl/2 



lhs. = C {N+l) ^L N ^' 2 I d[E\ 5(Ei) J] (E- + 

•> i=i i=i 

/ 1 iV+1 ~ \ 
lAW^exp^-— $> 4 2 j . (B.3) 

The power of the Vandermode determinant A^(x) = n i<7 -( 

Jacobian 

from the coordinate transformation. The constant C/ff+i),p arising from the group 
integration can be found in Mehta's book [17]. Now the integral over the (^-distribution 
can be performed. 

N 



lhs. = (N+ l)C (N+1)tP L N p I d[E] f[ \Eif(E% + 

J i=i 

lA^l^expU^^A . (B.4) 

We see that the resulting integral can be written as a GOE (GUE) average over N x N 
matrices. This is indicated by using Ei instead of E^ as integration variables. We go 
back to Cartesean coordinates U~ 1 EU — >■ H and find 

lhs. = (I + W<Z*M LNft (|det^(^ + z)~W) N 

°(A r ),/3 

= ^±^±M LNfj G N {z) . (B.5) 

This is the desired identity with L N/3 = Cn,p/(N + 1)CV +1 ^. Eq. (j4T|) is obtained for 

z = e + 2iv 2 k//3. 
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